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ABSTRACT 

Accurate and efficient numerical simulation of flows through turbomachinery 
blade rows depends on the topology of the computational grids. These grids must 
reflect the periodic nature of turbomachinery blade row geometries and conform 
to the blade shapes. Three types of grids can be generated that meet these min- 
imal requirements: (1) through-flow grids, (2) O-type grids, and (3) C-type 

grids. This paper presents a procedure which can be used to generate all three 
types of grids. The resulting grids are orthogonal and can be stretched to 
capture the essential physics of the flow. In addition, a discussion is also 
presented detailing the extension of the generation procedure to three- 
dimensional geometries. 
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BLADE GEOMETRY AND COORDINATES 

The development of the grid generation procedure begins by considering the 
electrostatic potential field generated by an infinite linear array of point 
charges in a two-dimensional space. The density of the charges is assumed to 
alternate between plus and minus 1. The mathematical expression for the com- 
plex electrostatic field is given by equation (1), where K(z - z 0 ) is the com- 
plex potential and i = \f-l. Eor a Cartesian space z = x 4- iy, z Q is the 
location of the zeroth charge and s is the physical distance between charges. 
For a blade row whose geometry is given on a blade-to-blade surface of revolu- 
tion, z = m + iO, where m is the meridional distance and 0 is the angular 
position; s then represents the angular distance between charges. 

t y 

O 


O +1 



— U ► 

t +1 

O -1 


x 


Array of point charges 


K(z - z Q ) 


00 

(-l) n In (z - z D - ins) 

n=-» 



( 1 ) 


BUde-to-blade surface of revolution, 
shoving m - 8 coordinates. 


AIMMii I I ttltfal i nil' liUiiiiiaii 


DEVELOPMENT OF EQUATION FOR POTENTIAL FIELD SURROUNDING THE BLADE 


The field equation for the array of point charges is expressed in closed 
form in equation (2) . Note that this expression is periodic in either the y- 
or 9-direction with period 2s. The potential field generated by distributing 
the fundamental solution (eq. (2)) over the surface of a blade is given by equa- 
tion (3), where Y(z 0 ) is the source density distribution on the blade surface. 
It is required that the real part of to (that is, F) be equal to 1 on the 
blade surface (eq. (4)). At large distances upstream or downstream of the 
cascade, £(z) is assumed to approach zero. This additional requirement is 
expressed by equation (5). Equation (3) evaluated on the blade surface forms a 
singular integral equation, equation (6), for the source density. Once y is 
known, the potential field surrounding the blade can be computed by direct 
integration of equation (3). 
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DEVELOPMENT OF EQUATIONS FOR SOURCE DENSITY DISTRIBUTION 


The singular integral equations for y can be solved by paneling methods 
similar to those used in solving potential flow problems in fluid mechanics. 

To employ these procedures, one first factors out of the integral equation its 
singular behavior. (For eq . (6) the factored form is expressed by eq . (7).) 

Next the blade surface L is divided into a series of segments. Over each of 
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polynomials in z Q . For the cases examined to date it was found that these 
terms could be approximated by a constant equal to their value at the midpoint 
of each segment. (This approximation assumes the length of the segment could 
be made smaller than the scale of the local geometric blade features.) Thus 
the singular part of equation (7) integrated over a segment is approximated by 
equation (8), while the regular part is approximated by equation (9). The 
auxiliary condition (eq. (6)) is approximated over a segment by equation (10). 
Upon introducing these approximations into equation (7) and restricting the 
value of z Q to the midpoint of each segment, a system of linear algebraic 
equations is obtained from which Y( z 0 ) can be determined. 
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BOUNDARIES OF COMPLEX ELECTROSTATIC FIELD SURROUNDING THE BLADE 


With Y known, the complex electrostatic field surrounding the blade sec- 
tions can be found. The real part E and the imaginary part p of the field 
form a periodic orthogonal body-fitted coordinate system. The contours 
E - Constant enclose the blade, while the curves n = Constant project from 
the blade (i.e., E ~ 1) to the periodic boundary (E = 0) . The curves extend- 
ing to upstream and downstream infinity are denoted q = 0, n = ±L, respec- 
tively. The locations of these bounding coordinate curves are found using a 
Newton-Raphson scheme with equation (3) to numerically generate the inverse 
mapping function. This procedure, however, because of its slow computational 
speed, was not used to construct the interior grid contours. 
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INTERIOR GRID 


The interior grid points are constructed from the solution of the inverse 
electrostatic problem in which (x,y) or (m,0) are specified as functions of £ 
or r) on the boundaries. The field equation for this is Laplace's equation in 
terms of (£>0) (eq. (11)). The solution of equation (11) which satisfies 
boundary conditions consistent with the coordinates of the bounding curvilinear 
curves yields the interior grid geometry. This solution can be obtained by 
either numerical or analytical procedures. 
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GRID CLUSTERING 


Grid clustering to capture the physics of the flow field can be introduced 
prior to or after the solution of the inverse problem is obtained. To insure 
orthogonality, the stretching functions used for clustering must be one- 
dimensional (i.e., eqs. (12)). The £ transformation can be quite arbitrary. 
For potential flow computations, a linear transformation is generally used. 

For viscous flows, one attempts to cluster grid points near the blade surface. 
An example of a transformation which can be used for this purpose is given by 
equation (13). The parameters m 0 ,mq control the degree of stretching in the 
transformation. The clustering of grid points in the n-direction requires 
special consideration to insure that grid point periodicity is maintained. A 
grid point located on the periodic boundary at n = D 0 has an image at 
T] = -n D . In order to maintain this property, the transformation in q must be 
an odd function of D over the interval -L to L. A simple transformation 
which exhibits this behavior is a polynomial in odd powers of n, an example 
being given by equation (14) . The parameters m 0 and mq in this transforma- 
tion are again used to control the degree of clustering. 
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GENERATION OF C-TYPE AND THROUGH-FLOW GRIDS 




Thus far the application of the grid generation procedure has been 
restricted to orthogonal O-type grids. To generate through-flow and C-type 
grids by the current procedure, the blade contours must be modified by append- 
ing slits of zero thickness to their surfaces. For C-type grids, one slit is 
used. Its origin is generally taken to be the trailing edge of the blade. For 
through-flow grids, two slits are appended, their origins being the leading and 
trailing edges of the blade. The shape of these appendages can be quite gen- 
eral. The generation of the grids associated with these modified blade pro- 
files proceeds as in the O-grid procedure outlined above. The generated grids 
are orthogonal and periodic. In the case of blunt blades, however, they 
exhibit a singular behavior at the slit attachment point. 
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GRIDS FOR CASCADES 


For cascades of nonzero stagger, the C-type or through-flow grids generated 
by the current procedure are discontinuous across the slit. This undesirable 
property can be corrected by a simple construction. First the two grid lines 
connecting the upper periodic boundary and the trailing edge and the lower 
periodic boundary and the trailing edge are found. Next the location and value 
of n corresponding to the periodic image of A on the lower boundary (i.e., 
A') and R' on the upper boundary (i.e., B) are determined. The spacing of 
grid points along the slit between EF is specified, which in turn determines 
the distribution of n along both sides of the slit. This in turn determines 
the distribution of n along the periodic boundaries AB and B'C' . The grid 
points along A'B' are required to be periodic images of the points along AB. 
This determines the distribution of n along A'B' and the grid geometry up 
to EB' . The construction of the grid downstream of BrC' proceeds in the 
same manner as outlined above, with the grid spacing along B’C' defining the 
clustering pattern. The resulting grids will remain orthogonal and periodic 
under this construction procedure. 


E 
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DEVELOPMENT OF THREE-DIMENSIONAL GRIDS 


Three-dimensional grids can also be developed by the current procedure. 
The geometry of a typical blade assembly as viewed in the meridional plane is 
shown in the accompanying figure. Let the surfaces of revolution describing 
the hub and shroud be denoted by r = f^(x), r = fg(x). A surface of revolu- 
tion bounded by these limits is given by equation (15) . Similarly, let m 
denote the percentage of distance measured from the leading edge of the blade 
along the surfaces of revolution (eq. (16))- On a surface of revolution the 
blade section geometry is given in terms of m and 6, where 6 is the 
angular location around the wheel. On the blade surface of revolution m 
and 0 can be expressed in terms of the coordinates of the current 

orthogonal system. The resulting coordinate system will be orthogonal 

on a £-plane and conform to the blade hub and shroud surface. 
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EXAMPLE O-TYPE GRID FOR TURBINE STATOR BLADE 


An example of an O-type grid generated by the current procedure is shown 
on the accompanying figure. The blade is a turbine stator with approximately 
90° of turning. No stretching was introduced in developing the grid. In 
regions of high surface curvature there is a high concentration of grid points, 
thus permitting accurate resolution of the local flow physics. Far removed 
from the blade, the concentration of grid points becomes sparse. This results 
in an economical distribution of grid points in regions of uniform flow. The 
grid as shown was generated to solve a potential flow problem. For viscous 
flow, a stretching of the £ contours would have to be introduced to capture 
the boundary-layer region. In addition, the r) lines would have to be 
clustered in the trailing-edge region to resolve the viscous wake. 


EXAMPLE THROUGH-FLOW GRID FOR A TURBINE ROTOR BLADE 




The next example is a through-flow grid for a high-reaction turbine rotor. 
This grid, as described earlier in the paper, was generated by appending two 
slits to the blade surfaces. The grid geometry is rectangular and periodic 
upstream and downstream of the cascade. Across the wake the grid geometry is 
seen to be discontinuous. The use of such a grid would require special care 
in transferring flow variables across the slit. 


EXAMPLE C-TYPE GRID FOR A TURBINE ROTOR BLADE 


The accompanying figure shows a C-type grid for the preceding turbine 
blade. Upstream of the blade the grid geometry is similar to an O-type grid, 
f while downstream it resembles a through-flow grid. For this grid the con- 

k struction procedure outlined earlier was used to insure continuous grid geom- 

fe etry across the slit. For potential flow analysis the grid as shown is quite 

= suitable. For viscous analysis a clustering of the grid lines near the blade 

S surface and trailing-edge region would have to be introduced to accurately 
resolve the flow physics. 
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EXAMPLE C-TYPE GRID FOR CASCADE OF NACA 0012 AIRFOILS 


The final example is a C-type grid for a cascade of NACA 0012 airfoils. 
This grid was generated by When-Huei Jou of Flow Research. The grid generation 
procedure as modified by him produces grids which are nearly orthogonal. The 
slight nonorthogonality is due to stretching functions used to maintain con- 
tinuous grid geometry across the slit and cascade periodicity. This grid was 
generated to solve a two-dimensional potential flow problem. It can also be 
used to solve a three-dimensional potential flow problem, provided the inviscid 
wake is convected downstream along the slit. For viscous flow calculation a 
clustering of grid points near the body and wake region would have to be 
introduced . 
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NACA 0012 airfoil. Stagger angle 30.00; pitch 0.794; 
cascade grid 128 x 16. 



